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Abstract 

^ . We present a detailed study of simple 'tree' models for off equilibrium dynam- 

. ics and aging in glassy systems. The simplest tree describes the landscape of 

0^ ■ a random energy model, whereas multifurcating trees occur in the solution of 



the Sherrington-Kirkpatrick model. An important ingredient taken from these 



. models is the exponential distribution of deep free-energies, which translate 

O ■ into a power-law distribution of the residence time within metastable 'valleys'. 



These power law distributions have infinite mean in the spin-glass phase and 
this leads to the aging phenomenon. To each level of the tree are associated 
an overlap and the exponent of the time distribution. We solve these mod- 
els for a finite (but arbitrary) number of levels and show that a two level 
tree accounts very well for many experimental observations (thermoremanent 
magnetisation, a.c susceptibility, second noise spectrum....). 

We introduce the idea that the deepest levels of the tree correspond to 
equilibrium dynamics whereas the upper levels correspond to aging. Temper- 
ature cycling experiments suggest that the borderline between the two is tem- 
perature dependent. The spin-glass transition corresponds to the temperature 
at which the uppermost level is put out of equilibrium but is subsequently 
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followed by a sequence of (dynamical) phase transitions corresponding to non 
equilibrium dynamics within deeper and deeper levels. We tentatively try to 
relate this 'tree' picture to the real space 'droplet' model, and speculate on 
how the final description of spin-glasses might look like. 
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Des journees entieres dans les arbres. (M. Duras) 



I. INTRODUCTION 

Aging experiments in glasses and spin-glasses |I|], 0, are now the focus of an intense 
theoretical |Q, 0, 0, 0], 0] and numerical |T0[, |n[] activity. Inspired by mean- field 



solutions of the spin-glass problem [0 , , a simple picture for the dynamics was proposed 
in 



14| , based on the idea that metastable states act as traps in the phase-space with broadly 
distributed trapping times. This picture naturally leads to aging (i.e. non stationary dynam- 
ics) and suggests phenomenological laws for the decay of the thermoremanent magnetisation 
(and of the a.c. susceptibility) which are in quite good agreement with experimental data. 

However, there are experimental features which are not compatible with the trap model 
proposed originally, which in fact corresponds to the phase space of the simplest kind - 
that of the 'Random Energy Model' (REM), for which Parisi's q{x) order parameter is 
simply q{x < a;) = and q{x > x) = 1 |]T5[. In that model, a metastable state is a 
single configuration; hence all dynamics is frozen if the system cannot "hop out". Stated 
differently, there is no 'bottom of the traps' dynamics (see Fig 1-a). 

This has various unsatisfactory consequences. For example, the equilibrium a.c. sus- 
ceptibility is zero for all frequencies in this model, in plain contradiction with experiments. 
This can easily be cured by allowing fast, small scale fluctuations to reduce the 'self-overlap' 
q{x > x) from 1 to a smaller value qEA (the Edwards- Anderson order parameter). This 



'dressed' REM behaviour was found recently in simple models for glasses [16|, where the 



'traps' can be very clearly identified, in particular in numerical simulations |T7 . 

More importantly, a REM landscape is insufficient to account for the subtle effects in- 
duced by small temperature cycling ||18|. Let us mention in particular the striking memory 
effect observed on the imaginary part of the a.c. susceptibility when the temperature is 
cycled as T — > T — AT — > T. The signal rises strongly when the temperature is first 
decreased, showing that new dynamical processes are restarted. However, when the tem- 
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perature is raised again, the signal recovers exactly the value it had before the period at 
T — AT, as if this period had not existed. A possible interpretation is that, in our 
language, there are 'traps within traps' (see Fig. 1-b): the intermediate period corresponds 
to non-equilibrium dynamics within a trap - with hops between 'supertraps' frozen out. 

The aim of this paper is thus to analyse in detail and generalize the REM-like trap model 
of ref. fl^ (corresponding to a 'one-step' replica symmetry breaking (RSB) scheme [l^) to 



a fully foliated tree structure ('full replica symmetry breaking' [T^). A large amount of pa- 
pers already studied various types of dynamics on trees [jl9|, directly inspired from Parisi's 



ultrametric construction. We however believe that our model is closer both to Parisi's con- 
struction and to reality (in particular because most of these studies analyse deterministic 
trees) . As in its REM version, aging appears most naturally within this framework and pro- 
vides an excellent fit of experimental data, which in turn allow to determine the structure of 
the q{x) characterizing the tree. We show that the data is already very well accounted for 
within a 'two-step' RSB approximation scheme. We discuss qualitatively T— jump experi- 
ments and 'noise second spectrum' within our model. We speculate on how this 'tree 
of states' could be interpreted in finite dimensional space, and rephrase the droplet model 
of Fisher and Huse |]2T]] in that context. We suggest that there is a spin-glass transition 



temperature associated with each length scale. We believe that our model, although still 
phenomenological, is of help to grasp the subtleties involved in the non-equilibrium dynamics 
of random or glassy systems. It also sheds light on the analytical results recently obtained 
for 'toy' microscopic spin-glass models exhibiting replica symmetry breaking 0,0. 

II. THE MODEL AND SOME ANALYTIC RESULTS 
A. The Single Layer Tree (REM) 



The simplest version of the model, already considered in [|1^, is defined as follows: the 



possible 'states' of the system live on the leaves of the one layer tree drawn in Fig. 1-a. 
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To each branch a is associated an energy barrier distributed as (minus) the energy in the 
REM, i.e.: 

P(e) = ;^exp[-f] (1) 

^9 ^9 

where Tg is the glass transition temperature. To each a thus corresponds an Arrhenius 

hopping time = roexp where Tq is a microscopic time scale. From Eq. (1), one finds 

that the Tq are (quenched) random variables distributed according to a power-law of index 
— r ■ 

^ rp . 

Pi-) = % > -o) (2) 

We shall call r an [x]-variable. For < a; < 1 (i.e. for T < Tg), the mean value of an 



[x]-variable is infinite: this will lead to the aging phenomenon ||14|| . 

The model is further defined by saying that hopping rate from state a to state j3 is -j^. 
This simple dynamics has the virtue of giving the correct Boltzmann equilibrium distribution 
: for finite A^, the equilibrium weight is: P^^ = Z^^ exp ^. The overlap between two states 
a,P is go = 0, while the self-overlap is taken to be qea < 1, independently of a. The 
'spin-spin' correlation function C(t^ + t,ty^) = {S{t + tw)S{tw)) within this model is thus 
naturally defined as: 

C{t^ + t,t^) = qEA^i{t,'tui) (3) 

where Ili{t,tw) is the probability that the process has not jumped between tw and t^, + 1 : 
this is the basic object we shall proceed to calculate. 

Note that from Eq. (3), C{tw,tw) can be smaller than 1. Strictly speaking, this is true 
in the limit tq = 0: it means that there is some dynamics taking place at short times, and 
refiecting the fact that a 'state' a is a ensemble of configurations mutually accessible within 
microscopic time scales. One could thus add branches connecting all these configurations 
to the 'node' a (see Fig 1-a); the 'trapping' times associated with this level of the tree are 
however such that their average is finite. A natural way is to continue Parisi's q{x) function 
for X > 1 and Qea < q < 1 |2^, such that the corresponding (r) = -j^tq are microscopic. 



This corresponds to the equilibrium part of dynamics (see p3[), while C{t, t^,) corresponds to 



the non-stationary (aging) part. This important idea will be generalized below: for a finite 
number N of branches, or an upper cut-off in the distribution (2), the distinction between 
equilibrium and non-equilibrium becomes itself time-dependent. 

The system is assumed to start in a randomly chosen initial metastable state. Rather 
than work at fixed tyj, it is convenient to work in Laplace transform, and define ni(t, E) = 

Edtyj exp{—Etyj) Ili{t,tw). This is because, conditional on the value Tq,, the amount 
of time spent in a trap is exponentially distributed - using some standard properties of 
exponential random variables then makes calculations relatively straightforward. 

We define p{E,Ta) to be the probability that at a random exponential time of rate 
E the system is found in trap a with hopping rate — . To be in trap a at time t^j, one 
possibility is that the system starts in trap a and time occurs before the system escapes 
from that trap. Clearly one starts in trap a with probability the probability that iw 
occurs before escaping trap a is then given by property (3) of exponential random times as 
E/(E + —). However, to start with, the system may be in any of the traps 3 1 < < N 
but escape from it before time t^. Conditional on this having happened, the memoryless 
property of the exponential time (property (A-1)) means that the probability of being in 
trap a at time is just p(£', r^), that is the game starts anew. Putting all of this together 
allows us to write down the renewal equation 

hence 



In the case where N is large we may use the approximation 

p{E,Ta)= ^""t^ (6) 

the angled brackets indicating averaging over the disorder. In the regime where Etq -C 1 
(equivalent to ty, ^ tq) we obtain 

(-^1^) ~ xr^E^cix) (7) 
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where 



TT 

c(x) = r(a;)r(l — x) = — ; — - 

sin(7rx) 

We now wish to calculate the probability that after an initial waiting time the system 
has not decayed from the state it was in at time after a subsequent time t has elapsed 
after tw Keeping the waiting time exponential, we find 



N 



Il,it,E) = {J2piE,T^)e 



(8) 



Now, making a judicious change of variables and noting that exp(— i^w) = /q°° Edv exp{—Ev) 
6{v — w), one finally finds: 



Sm TTX 



TT 



x—l„—x 



du{l - uY u 



(9) 



This last form shows very clearly that the correlation function only depends, within this 



model, on the rescaled time i.e t over the 'age' of the system t + 1^: as argued in [jT4 



only traps such that r ~ 0{t.u}) have an appreciable probability to be observed after time 
tyj. From the same argument, the average energy decreases, in this model, as £{t^) — £{t = 
0)~-Tlog(^). 

Two asymptotic regimes are of interest. If t -C one obtains 



ni(t,t. 

In the regime t^, we obtain 



1 - 



Sm TTX 



ni(t,t. 



vr 1 -x) \t + t 



sin(7rx) / t 



l-x 



(10 - a) 



Tlx \t + t 



(10-6) 



C{ty^ + t,tw) has thus the 'weak-ergodicity breaking' property [|T4 



but 



lim C{tyj + t,t^) = for all finite 

t — too 



lim C{tu] + t, tw) = qEA for all finite t 

^oo 



;ii-a) 



:il-6) 



meaning that 'true' ergodicity breaking only sets in after infinite waiting times |^ 



The critical case x = 1 is of special interest, since it leads to aging but with an 'anomalous' 
scaling (i.e. not as 7^), we find: 



hn 1 — exp ( - 

ni(t,t^) =To r du- "° ( (12 - a) 

Jo {logu){tw + t - UTq) 



tw+t—UTQ 
TO 



which leads to 

ni(t,tj = 1 - 1^ for log(^) « log(^) (12 - b) 

and 

ni {t, tn,) = tt\^ for ^ > (12 - c) 

The above results pertain to the 'aging' case x < 1. For completeness, and also because 
they may be relevant to describe the equilibrium dynamics at short times (or for T > Tg, 
see [0), let us give the corresponding results for x > 1 in the limit ^ tq: 

Ui{t, t^) = l ^(-) for t < To (13 - a) 

X-1 To 

ni(t,t^) = (x-l)r(x)(^r-i forro<t<t^ (13-6) 

V 



Mt,U = i^- l?r{x)i^^-^) for t » (13 - c) 



and 

Note the major difference between these last results and those given in Eq. (10): for tq — > 0, 
C{ty^ + t, tw) is identically zero for a; > 1, but remains non zero for x < 1. In this sense, the 
crossing of x = 1 corresponds to a true dynamical freezing transition. 

B. The Multi-Layer Tree (RSB) 

The above single layer tree model can be generalized, following the construction of 
Mezard, Parisi and Virasoro for the tree of states (see [|12|, and also p6[). The hopping 



rate from state a to state j3 is set by a random escape time Ta^j^ the statistics of which 
depends on the overlap g^,/? between states. In Parisi's tree, the free-energies associated to a 
certain level are distributed according to Eq. (1), but with a level-dependent Tg, and hence 
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with the overlap associated with that level. Ta,p is thus an [x]-variable, with a gr— dependent 
X which is precisely the inverse of Parisi's order parameter function q{x). In the case of the 
REM {q[x <x)=Q and q{x > x) = Qea) we recover the single-layer tree since all the r's 
have the same statistics. 

The procedure we follow to solve this model for a finite number M of layers (finite RSB) 
is a generalisation of the one presented above. We introduce the probabilities Ilj{t,tu;), 
j — 1, M that the process has never jumped beyond the j^^ layer of the tree between t-w 
and tw + t. j = M corresponds to the deepest level of the tree and j — I the upper level 
of the tree - see Fig 1-b - and hence Hi is the same as the one introduced for the single 
level tree. Introducing the j*^ level overlap qj (with go = and qM — qEA), the spin-spin 
correlation function is clearly: 

M 

c{t^ + 1, t^) = ij [nj(*, tw) - n,+i(i, t^)] (14) 

with the convention that no(t,tu,) = 1 and YlM+iii-iiw) = (all processes inside a state 
have indeed taken place on the time scales ^ Tq - see the discussion after Eq. (3).) The 
precise construction of the multi- layer tree is a follows, the {j + l)th layer tree is constructed 
by adding Nj^i branches at the end of each branch of the j layer tree - the new states 
are now at the ends of these new branches. Associated with each new branch a is a new 
hopping rate — - — with which the process hops back to level j and then falls back into 
the (j H- l)th layer. For simplicity, we define this hopping to be independent of all other 
hopping on the tree, i.e. the hopping dynamics is in parallel. We therefore have transitions 
from states characterised by inverse hopping rates {ti, T2, • • • , Tk-i, Tk, • ■ ■ , tj-i, tj} to states 
{ti, T2, ■ ■ ■ ) Tfe, ■ ■ ■ , 'Tj} ^-t rate (the primes denoting new rates chosen within 
the tree structure) and all these transitions are in parallel. The tj are quenched random 
variables chosen according to the law 

Pir,) - ^ (T, > TO,) (20 

T- 

where < Xi < X2 ■ ■ -a^j-i < Xj < 1. These times Tj corresponds to (free) energy barriers 
ej distributed as in Eq. (1), but with a level dependent freezing temperature, Tgj. Hence, 



again, the dynamics of this model reproduces, at long times and for a finite number of 
branches, the Boltzmann equilibrium: the total weight of a branch at level j is given by 
exp ^. 

For technical reasons we shall also assume that microscopic time scales associated with 
levels lower down the tree are much smaller than those associated with microscopic time 
scales further up the tree - this is also a reasonable hypothesis from the physical point of 
view. Therefore we shall impose Toj+i <S Toj. 

Once again we work with an exponential random time tyj of rate E. We define 
Pj{E; Ti, • • • , Tj) to be the probability of finding the system in a set of traps characterised by 
{ti, ■ ■ ■ ,Tj} at time tw We calculate the pj inductively as follows 

Pj+i{E; Ti, • • • , Tj+i) = pj{E; Ti, • • • , Tj) ■ p{Ej+^,Tj+i) (15) 

where Ej ^ E + Ei^i ^/n for j > 2, Ei ^ E and p is the result for the single layer tree. 
This comes from an application of Baycs' theorem and using properties (A-2) and (A-3) of 
exponential times. Consequently, Pj{E; ti, • • • , Tj) = Y[i=iP{Ei, Ti). The fact that we have 
chosen tqj <S To(j_i), and in addition chosen M to be reasonably small, means that we are 
assured that EjTqj <^ 1 and hence we may use the small E asymptotic result for the p{Ej, Tj) 
in the expression above, yielding 



p.(^;n.-.^)-n ^,.^;^,.^ (16) 

We now use 

fi,(t,E)= pAE;n,---,T,)expi-tj2-) (17) 

{ri-r,} i=l 

and the fact that the A^j — > oo and hence all combinations of {ti, ■ ■ ■ ,Tj} exist, allows us to 
write 

J 1 \ J i 

-^H-m^B^. (18) 

i=l J i=l C\Xi)l^i 

One again attempts an integration by substitution, starting with Tj and proceeding down 
to Ti. On inverting the Laplace transform one obtains 

nLB(.<,i-x.) '''' 
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where 9 is the Heaviside step function. The above formula has an immediate probabihstic 
interpretation. If Ui are independent random variables from Beta distributions of indices 
1 — Xi and Xi respectively, then 

Uj{t,tJ = P{flU,>-^) (20) 

A useful rearrangement of equation (19) for asymptotic analyis is 

Jo' ni=i duiu-'^'ii - Ui)-*-'eini=i Ui - s)b_^{i - x.+o 

Ilj + l{t,tu)) — Ilj{t,t.w) 



where s = tvt- and 



ng5(a;„l-a;,) 

(21) 



(x,y)= r duw'-^il-uy-^ (22) 
Jo 



is the incomplete Beta function (clearly Bi{x,y) = B{x,y)). The above rearrangement also 
makes clear that the obvious inequality Ilj+i{t,tw) < Ilj{t,tu,) is satisfied. 

Let us focus on the short time and long time asymptotics implied by the above analytical 
forms. For t ^t^, we find that 

r(f^ff)-n in a \ ^^"(^^^) ~ ^^'^ i ^ \ I 

7r(l - xm) jj[ B[xj, 1 - Xj) \t + ty,J 
+0((l)i-^M-i) (22) 

This form shows that the initial decay is, as expected, dominated by the deepest level of the 

tree. The possibility of 'jumping' to higher levels during the waiting period however acts to 

1 

reduce the effective waiting time by a factor n?li^ Bixj,XM-xj) 



^J=l B(xj,l-Xj) 

large scale jumps indeed completely restarts the small scale dynamics. 

In the limit t t^j, on the other hand, the asymptotic decay is given by: 

sin(7ra;i) f Y\ r^u^y^Mx^^ 



^ < 1. Each of these 



C(t. + t,U = ?i— hr^ TXT +0{{^r^^') (23) 

TTXl \t-\-ty,J t 

Note that, again, C{ty, + t, tj) is only a function of the ratio j-. This would not be true if 
the number of branches Nj is finite, as we shall discuss now. 
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C. Finite number of branches and interrupted aging. 



Let us now discuss the case where the number of 'states' is finite, starting with the 
single layer tree model. As can be seen from Eq. (5), or directly from the argument giving 
the order of magnitude of the largest r drawn from distribution (1), the results established 
in section II. 1 rely on the inequality: 

t + < ToN^ = terg (24) 

In the opposite case, the dynamics ceases to depend on t, t^, and stationary dynamics resume 
and aging is 'interrupted': tgrg is the ergodic time of the problem. Hi asymptotically behaves 
as in Eqs. (10), but with the variable replaced by t^. The scaling of C(ty; + t,t^) as 
a function of -f- is thus expected to degrade progressively as approach terg(A^)- This 
was proposed as an explanation for the small (but systematic) deviation from a perfect 
scaling observed experimentally [1^. Let us note that very similar results would of course 



hold if N is infinite but the power-law distribution (2) happens to be truncated beyond a 
certain terg- 

The case of a multilayer tree is interesting, but rather complex - we shall thus only give 
qualitative arguments. Let us assume that the ergodic time at level j is now terg"'- Is is 
reasonable to assume that the hierarchy of the Xj (slower dynamics as one gets higher up 
the tree) is maintained for the ergodic times, i.e.: terg'' ^ ''^erg'''''^- For a given waiting time 
tu,, a specific level of the tree is naturally selected through the condition 



''erg ''erg \ 

meaning that all the levels 'below' j{tw) are equilibrated, while those above are still aging. 
In this case, even if < 1, the corresponding transitions are equilibrated - reducing 

further the effective value of QeAi renormalized by small scales dynamics. The interesting 
point, however, is that the effective 'ergodic' time for which deviations from a scaling 
first become visible is now tgrg"'^*™'* which becomes, through Eq. (24), dependent on t^ itself. 
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[It might well be that the ergodic time determined experimentally in is only the one 
corresponding to the particular level of the tree which happens to be probed on the 1 — 10000 
seconds time scales]. 

Another consequence of these nested equilibrium time scales is the following: if a per- 
turbation is applied to the system at time t^, the response at time + t will depend on 
how much the system has been able to reequilibrate during time t. The response is thus 
expected to depend on the particular level j defined as in Eq. (25), but with t replacing 
tw Hence, for t <^ t^, one expects that C(tw + t,tw) = qj{t), leading to violations of the 

scaling at short times too. Said differently, the effective 'microscopic' time scale tq is 
replaced by terg-'^^'^^'^^ '^^^ existence of different time domains in the plane tw,t was in 



fact suggested in refs. 0, 

III. RELATION WITH MEASURABLE QUANTITIES. 

A. Constant temperature data 

Although noise measurements can be and have been performed on spin-glasses to access 
directly the correlation function computed above, the response function - measured through 
the thermoremanent magnetisation (TRM) or the a.c. susceptibility - is of much easier 
access. Most of the available data thus gives information on the response function R{t, t') = 
^dHfJ)' (^^^ below). For equilibrium dynamics, the response and correlation functions are 
related through the Fluctuation-Dissipation theorem (FDT). For non-stationary dynamics, 
an extended form of this theorem has recently been proposed HH, [0, 0, This relation 
can be written as: 

T R{t,t') = d{t - t')X[C{t,t')f^^^ (26 - a) 

where X{.) is a certain function which can be computed within simple models P|, 0, [§. 
For the SK model near Tg and the toy-model of X is the inverse of the function q{x), 
whereas for the spherical p— spin model, X is a constant < 1 (depending on temperature). 
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(Equilibrium dynamics corresponds to X = 1 and one recovers the usual FDT.) 

Within the present context, the response function can be computed by assuming that to 
each trap i is associated a certain magnetisation rrii with probability p{m), which we shall 
take to be even p{—m) = p{m). The hopping rate from state i to state j is modified in the 
presence of a field H by a factor exp — ^i'^'^^-^-O'^-j) ^j^i^]^ ensures the correct equilibrium 



weighting for all values of C. It is then easy to establish that 

T R{t,t') = e{t-t' 



dt' ' dt 



(26 - h) 



If C{t, f) = C{t — t'), the usual FDT is recovered. If - as is the case here - C(t,t') = C{jj 



t' n 



then one finds that Eq. (26-a) holds with X = (1 - C) + Cy = (1 - C) + c^%)- ^ote 
that in the early regime corresponding to equilibrium dynamics {t = t' + 0(ro)), X = 1 
independently of C- 

Now, the TRM M{t^ + t,t^) is defined as: 

M{tu, + t,tu,) = dt'R{t^ + t,t') H (27) 
Jo 

where H is the applied field between and t^. Defining the function Y{q) through X(g) = 
^^^^ and using Eq. (26), we obtain that the experimentally measured TRM and C{tw + t, t^) 
are related through 0: 

M{t^ + t,t^) =Y[Cit^ + t,t^)]h (28) 

where h = The simplest assumption of a REM with q{x < x) = 0, q{x > x) = qea and 
X = const, leads to: 

M{t^ + t, t^) = QEAXhlliit, t^) (29) 

Experimentally, the decay of M{ty^ + t,tw) at short and large times is indeed very well 
described by the forms given in Eqs. (10-a, 10-b). In particular, the scaling as a function 
of the reduced time is reasonably well obeyed (but see ||14[-b] and section II. 3 above). In 
fact, the experimental value of x is in the range 0.6 — 0.9 if extracted using the short-time 
expansion (10-a), and in the range 0.05 — 0.2 if extracted using the long-time expansion (10- 
b) fl^, where it was called 7]. This suggests that a more complicated RSB scheme is needed, 
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generating a tree with at least two levels - allowing for the existence of two independant 
exponents Xi and X2. In the spirit of Parisi's step by step construction of the SK solution, 
we shall thus assume that q{x) is given by: 

q{x < Xi) = q{xi < x < X2) = q q{x > X2) = qsA (30) 

In order to obtain M{tw + t, we first assume that X = const., or ( = 0. This assumption 
is natural within our picture since qj corresponds to the fraction of 'frozen' spins at level 
j — 1: these spins thus retain their magnetisation when the process 'jumps' at level j — 1 
leading to: 

M{t^ + t, t^) (X [{qEA - g)n2(t, K) + qlllit, t^)] (31) 

This form is only valid after the short time (x > 1) relaxation has taken place; just before 
the field is cut, the field-cooled magnetisation is by definition M{t^,tw) := xfcH, where 
Xfc is roughly temperature independent. We choose units such that M{t~,tw) = 1. The 
four fitting parameters are thus Xi, X2, for the shape of the TRM decay and for the 
?/— axis scale. We show in Fig 2-a and 2-b the best fits obtained for both studied samples, 
Ag:Mn2.6% (AgMn) at T = (0.87 Tg) and CdCr1.7Ino.3S4 (Crin) at T = lO/sT (0.6T,,), 
plotted as a function of the rescaled variable s = j^j-. As can be seen, these fits are excellent 



on the whole time regime, which extends, for CrIn, from 10 sec. to 2.3 x 10^ sec |^9|. The 
value of the fitting parameters for different temperatures and samples are given in Table 1 
and 2. Assuming that X is not critical at Tg, we can extract values for the exponents 
defined as qEA — Y and q ~ (^^" Y' : j3 ~ 0.65 and ~ 1.0 (both samples give 

the same values for these exponents, within 10%.) Our value of (3 is compatible with other 



determinations (/? ~ 0.5 — 1) |^ 



The fits shown in Fig. 2-a, 2-b implicitely assume that 0:2 < 1. We have also tried the 
case X2 = 1 - the free parameter now being tq. The resulting fit is of much poorer quality; 
however, as will be discussed in section IV (Fig 6 below), the initial part of the TRM decay 
for different is rather well rescaled when plotted as a function of as suggested by 

Eq. (12-b). We interpret this as a sign that q{x) is in fact continuous, at least near x = 1. 
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We have also tried to fit the data using different values of (. The quality of the fit remains 
excellent if C < | (although the values of qEA,(l are somewhat changed), and deteriorates for 
larger values of (. One can also use the form of X{C) suggested by the models studied in 0, 

0, 0, i.e., the inverse of the function g(x) defined by Eq. (30). Defining a characteristic 
time t* such that C{tw + t*, t^) = q, the TRM takes the following form (for t = 0+): 

Mit^ + t,t^) = [xi{Cit^ + t,t^)-q)+X2q]h (t < t*) (32 - a) 

M{t^ + t,t^) =X2C{t^ + t,t^)h it>t*) (32-6) 

1. e., the decay rate is discontinuous at t = t*. We have plotted in Fig 3 the best fit obtained 
using Eqs (32): the initial and final part of the curve can be accounted for, but not the 
intermediate region around t = t*. This suggests that either the function X is (in our case) 
not connected simply to q{x) and is a constant for the range of correlations probed, or else 
that a continuous RSB scheme is needed (to remove any discontinuity in decay rate). More 
work on this point is certainly needed, in particular to extend our formulae to a continuous 
tree. 

Let us now turn to the a.c. susceptibility measurements. The object which is measured is 
now xi^i ^) = /-oo ^(^^ Gyip{iij{t' — t))dt' . Let us first assume a REM-like dynamics. Then 
using Eqs. (26), and in the limit t ^ uj^^ (actually needed for an accurate measurement of 
x{tj,t), we find that: 

^(co, t) = r du (exp[-z(a;t)M] - 1) ^^^(^^) ^-^ (33 _ a) 

T Jo vr 

(the = value of x' has been removed). Defining uJ = ut and changing variables, 

x(o., t) ^ xiuJ) = x^^^r(l - x)e-^uJ--' (33 - b) 

in agreement with the arguments given in 11. Eq. (33-b) shows that: 



(i) the relevant variable is again the rescaled time ujt: this is very well obeyed by exper- 
imental data. 



16 



(ii) There is a 'constant loss angle' S = |(1 — between the aging part of x' and x", as 
expected from Kramers-Kronig relations. 

(iii) For a multilayer tree, the same formula (33-b) will hold in the limit uJ ^ 1, with x 
replaced by xm- In other words, the a.c. susceptibility measurements are only sensitive to 
the deepest level of the tree. Actually, Eq. (33-b) for x" is very well obeyed experimentally 
- although a 'non aging' contribution Xeq. must be added: see Fig. 4. Indeed, in the limit 
uj — > oo, the short time dynamics corresponding to intra-state fluctuations {x > 1) cannot 
be neglected in Eq. (33). 

As mentionned above, magnetic noise measurements have also been performed in order 
to test FDT in spin-glasses The 'usual' FDT would predict that the noise spectrum 
S{u!,t) is related to x"(u;,t) through: 

X'\uj,t) = - ^ (34) 

with X = 1. The non-equihbrium theory proposed in [|], |^, claims that X = xm < 1 
in the limit ujt ^ 1, while our model suggests (in the same limit) X = 1 — ^ - see also 
||32|| . Experimentally, the noise measurements were not calibrated in an absolute way. The 
proportionality relation (34) was indeed confirmed but no value of X could be determined. 
In view of the present discussion, it seems to us that this would be a very important point 
to check: the value of xm extracted from the decay of x" the TRM should be the same 
as the proportionality constant extracted from Eq. (34). Since xm is in the range 0.6 — 0.9, 
a significant effect could show up by comparing carefully equilibrium and non-equilibrium 
regimes. 

Let us finally note that for a; > 1, x"(u;,t) is given by: 

X"(w, t) oc for u < t-^ (35 - a) 

X"{uj, t) oc {uoTof-^ for < < t^^ (35 - h) 

(the large frequency region uj ^ t^^ requires a detailed description of p(r) for r ^ tq). 
Hence for x ~ 1 (i.e slightly above the spin glass transition, or equilibrium dynamics in the 
spin-glass phase), Eqs (33-35) gives S{ijj) oc ^ i.e. the model naturally leads to 1// noise. 
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B. Temperature cycling 



As mentioned in the introduction, more sophisticated experimental protocoles have been 
investigated [0, |3^, where the temperature or magnetic field is changed during TRM 
or x" measurements. A complete analysis of these experiments is beyond the aim of the 
present paper and will be given in separate publications |Q. We shall restrict here to a 
qualitative discussion of the results, focusing on the remarkable memory effect mentionned 
in the introduction, and illustrated in Fig. 5 (from |p5|| ). 

Let us thus discuss the effect of a small temperature change on x" within a REM land- 
scape. As discussed in the previous section, this REM description is sufficient to account 
for the behaviour of x"{^,t) in the limit ut ^ 1 (provided an equilibrium contribution is 
added). In the REM, the (free-) energy landscape does not evolve with temperature and x 
is simply given by A small change of temperature T — > T' = T — AT simply changes 
the trapping times as: r — >• r' = ro(— )^ with p = ^ > 1. Suppose that the waiting time 



before the small temperature jump is t^i - see Fig 5. As emphasized in [|1J] and above, only 
traps with r ~ t^i have a significant probability to be observed. Hence, if the time spent at 
temperature T' is such that tw2 ^ 'T'oi^)^, it is quite obvious that no jump will take place 
during this intermediate period. The 'memory effect' x"('^5^u'i + + 0) = x"{^:'twi — 0) 
would then be trivial. On the other hand, one can show that the change of x" at t^i is given 
by: 

x"{uj,t^, + 0) qEAiT')X{T') , 
X"{uj, tyji - 0) qEA[T)X{T) 

Typically (see Fig 5), p = 1.2, x = .75, while the variation of qea is roughly a factor 1.2. 

Hence, provided tq < 10~^ sec, x" should first decrease and then remain constant when the 

temperature jumps down, at variance with the experimental results. In physical terms, all 

dynamical processes should be slowed down by the temperature jump - hence the decrease 

of x" and the memory effect. The simultaneous observation of a strong increase of x" and 



dx" 

at 



['rejuvenation') and of the memory effect is non trivial In |T8[, it was suggested 



that the origin of this combined effect comes from the hierarchical nature of the phase- 



space, with valleys bifurcating into subvalleys at all temperatures below Tg in a continuous 
sequence of 'micro phase transitions'. Even a small temperature change can thus be thought 
as a 'quench' from 'high' temperature, setting back the age of the system to zero - hence 
the increase of x"- 

We want to rephrase this interpretation within the framework developed here. As men- 
tionned above, the construction of Parisi's tree can be extended 'within states' to x > 1 and 
Q > Qea, simply reflecting the fact that even within each state there are long-lived configu- 
rations and non exponential equilibrium dynamics, even above Tg. However, as mentionned 
previously, only the part of the tree corresponding to x < 1 corresponds to non- stationary 
dynamics, while the part x > 1 is equilibrated within microscopic time scales. Tg is the 
temperature at which the smallest x first reaches the critical value 1, but one can expect a 
whole sequence of phase transitions as the temperature is decreased, corresponding to the 
crossing of x = 1 for the successive levels of the tree, much as in the 'generalized' Random 
Energy Model (As discussed below, this could correspond in real space to the progres- 
sive freezing of the dynamics over smaller and smaller length scales.) Hence if at a certain 
temperature xm < 1 but xm+i > 1, only the transitions at the M**^ level of the tree will 
contribute to the relaxation of x"- But if upon lowering the temperature xm+i becomes 
less than 1, then the processes corresponding to this level of the tree will start contributing, 
while those corresponding to xm are to a certain extent frozen - allowing for the memory 
effect. The small temperature jump is furthermore tantamount to a quench for the M + 1 
level since when xm+i > 1, all states are more or less equivalent, corresponding to a random 
initial condition at t = t^i + 0. 

Another related scenario is also possible: as discussed in section II. 3, the finiteness of the 
ergodic time terg"* at level j leads to 'interrupted aging' 
levels of the tree are in this case not those such that xa/ < 1 < xm+i but, as discussed 
in section II. 3, those such that tcrg^^ (T) ^ ^ tcrg^^^^iT), even if xm+i < 1- Aging 
processes will be restarted upon cooling if terg*^^^(7' — AT) ^ t^. A quantitative analysis 
of the experimental data along these lines is currently underway. Experiments investigating 
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, after a time terg"'- The crucial 



the role of the magnetic field on aging and their interpretation are reported in 



P7| 



C. Noise second spectrum 

Another interesting set of experiments which can be discussed within our model is the 
noise 'second spectrum' measurements of Weissmann et al. Very briefiy, the point is 
that on sufficiently small samples, one observes j noise, but with an amplitude which is 
itself 'noisy', i.e. randomly changing with time (albeit on rather long time scales.) The 
natural interpretation - discussed in different terms in pO| - is that the noise primarily 



comes from the near equilibrium levels x ~ 1, but that higher level jumps will slightly 
change the amplitude of this noise - since the systems are small enough, the averaging over 
different 'branches' is not performed. For a multilayer tree, the prediction is obviously that 
the noise spectrum behaves as uj^m-'2 ^ u^'^, but with a noisy amplitude with a correlation 
function decaying as t~^^, as observed experimentally (where this exponent is called (3). A 
more quantitative analysis of the experimental results along these lines would however be 
desirable, but we note that the value of P = Xi determined for CuMn in is very close to 
the one quoted in Table 1,2 for xi, and evolves similarly with temperature. 

IV. DISCUSSION. TREES AND DROPLETS 



Fisher and Huse |^T|], see also were the first to stress the importance of understanding 



the nature of the spin-glass excitations in finite dimensional spaces. They argued that these 
excitations should be of the form of compact 'droplets', with a surface much smaller than 
their volume. These droplets were defined with respect to the (unique up to a global spin- 
fiip) ground state as the 'spin-fiip' excitation of lowest possible energy within a region of size 
L] the energy of such a droplet is found to be of order {6 ~ 0.2 in = 3), and the energy 
barrier for such a droplet to be activated grows as L"^ . Apart from the presence of these 
ever-growing energy barriers, the description of Fisher and Huse of a spin-glass is that of a 
'disguised ferromagnet'. In particular, exactly as in a ferromagnet quenched in zero field, 
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a spin-glass quenched below Tg would approach equilibrium by growing domains of the two 
'pure phases' assumed to be present, thereby eliminating the extra energy (L^) associated 
with the domain walls. The system thus coarsens with time; after time tw the typical size 
of the domains is R{tw) = log'''(^) - instead of R{tw) = tU"^ in a ferromagnet (with non 
conserved dynamics, see |^). Before discussing the form of the relaxation proposed in ||21[ 



inspired from this picture, it is useful to recall how C(t^ + t, tw) would look like for a simple 
ferromagnet. An exactly soluble case is the 0^ theory in the limit n — > oo, for which one 
finds 



Cn=oo{t^ + t, tw) - i (2il + t)2") ^^^^ 



(d is the dimension of space). Mazenko's approximate theory for n = 1 [Q, |39| leads to 

4 



Cn=i{tw + t-w) — tan 



-1 



TT 



2 Cn=oo(tw '^tyt'w) -. foo\ 



For t ^ tu,, expressions (37,38) simplify to [^)'^/^ = {^^^^)^ with A = |, whereas more 
elaborate theories show that A is non trivial in general [^. Fisher and Huse argued that 
the appearance of the ratio of two length scales, R(t) and R(tw), should be general and 
proposed for the spin-glass dynamics to write: 

CpHit^ + t, tw) oc i?(t)-^S(^^) (39) 

n{twj 

with S(a; — > 0) = const, and T,{x — > oo) = x^^'^. (Again, we identify C(tw + t,tw) and 
M{tw + t,tw) in the small H limit). Eq. (39) means that: 

log(^) 

a) The scaling variable for aging should be - — tSt rather than 7^. 

b) Even in the limit where t^ = 00, CpHitw + t, t^) decays to zero at large times if 6' > 0, 
at variance with Eq. (11-b) and the results of P, |0, 0. 

Before comparing with the experimental data, let us remark that, from Eqs. (37-38), 
Citw ~l~ tw ) is a function of over the whole range of time scales. This scaling variable is 
thus more general than the asymptotic one, i.e. -^^^ This suggests that even in the case 
of logarithmic domain growth, a 7^ scaling should hold for coarsening. We have confirmed 
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this for the Id Random Field Ising Model where R{t) oc log'^t (see e.g. [0], ^^). As 
emphasized in [|], Eqs. (37-38) show that aging is not specific of spin-glasses, but occur as 
soon as the equilibration time of the system is infinite (or much larger than t,tw)- We have 
plotted in Fig. 6 Cn(tw + t,tw) as given by Eqs. (37-38) as a function of j-^- The major 
difference with the spin-glass data is the singular behaviour of the latter for short times (see 
Figs. 2-a,b). 

We have tried to test the predictions contained in Eq.(39). For example, we have looked 
for the best values of 6 and tq to rescale TRMs at a fixed temperature and different t^, 
for both AgMn and Crin. The data points unambiguously towards ^ = 0; data collapse is 
quite good at short times but deteriorates as soon as t > : see Fig. 7. The other test is 
to compare different temperatures with the same tw We fixed ^ = 0, and chose to(T) and 
the T-dependent proportionality constant in Eq. (39) (which is essentially qEA{T)) to best 
rescale the late part of the curves. Although this late part can indeed be satisfactorily fitted 
by log~^(;^^^^) with ^ ~ 1.2, the early part of the data does not scale (Fig. 8). Hence we 
believe that Eq. (39), even with = 0, is inadequate to describe the data consistently. 

Koper and Hilhorst ||3^ proposed a similar picture, although they assumed that R(t) 
grows as a power-law rather than logarithmically. Also, the relaxation time associated to 
a domain of size R is taken to be oc i?^. It would be too long to discuss in detail this theory 
here; let us simply mention that it does lead a relaxation of x"{^,t) behaving as uj~^t~P^. 
Consistency with the ut (or ^) scaling observed experimentally leads to p^; ^ 1 - which 
is in fact natural since it means that the growth time and relaxation time of droplets are 
comparable However, the resulting relaxation of x" as is much too fast (compare 
with Eq. (34), with x ~ .75). 

The basic remark of Fisher, Huse, Koper and Hilhorst that the flipping spins are clustered 
somewhere in space and that the time scales should grow with the size of these clusters 
seems however unavoidable. How can this be reconciled with replica symmetry breaking ? 
Of course, a proper replica (or dynamical) theory in finite dimension is needed to answer 
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completely this question. Such a theory is not yet available for spin glasses, but has been 
worked out for the simpler problems of manifolds in random media (i.e. polymers, 



surfaces, vortex lattices, etc.). One basically finds that q{x) becomes L (scale) dependent, 
with a characteristic value of x{L) varying as L~^: small scales correspond to large x. 

From Eq. (1) and the interpretation of x{L), one thus sees that the energy distribution 
for the excitations ('droplets') of scale L has the form: 

1 Ifl 

p(/,L) = -^exp(--^) (40) 

showing, as postulated by Fisher and Huse, that the energy scale grows as L^. At this stage, 
an important difference with Fisher and Huse is the exponential form of the distribution. 
Interestingly, Eq. (40) means, within our interpretation, that there is a spin-glass transition 
temperature Tg{L) associated with each scale, defined through x{L,Tg{L)) = 1. The infinite 



sequence of micro-phase transitions suggested in 0, [|T8[ thus corresponds to a progressive 
'weak ergodicity breaking' (in the sense that x{L) crosses 1, see II. 1) of smaller and smaller 
length scales (faster and faster degrees of freedom). 

Quite naturally, one can associate an L-dependent overlap q{L) between configurations 
obtained by flipping a droplet as 1 — g(L) oc (j)'^-^ , where df is the fractal dimension of the 
droplets {df is equal to d for Fisher and Huse) and ^ is a correlation length. Hence, the 
picture we propose is in fact very close the original droplet model, except that: 

a) The 'droplets' are only required to be metastable, and not 'lowest' excitations, 

b) These droplets may a priori be non compact {df < d), and 

c) The time scale associated with droplets of size L is not peaked around some t{L), but 
rather power-law distributed with a parameter x{L) which becomes smaller and smaller at 
large sizes. 

d) The Fisher-Huse time scale Tq exp(L^) could be interpreted as the ergodic time terg(-^) 
associated to scale L: one indeed should expect that Eq. (40) ceases to be valid for |/| >>> 

e) Aging will be totally interrupted when the 'terminal' ergodic time scale terg(0 associ- 

23 



ated to ^ is reached. 

This picture suggests that the final description of the experiments should involve a con- 
tinuous tree, rather than the 2-level approximation that we have chosen. Further work on 
this aspect and on the related problem of the Fluctuation Dissipation theorem is certainly 
needed; we however hope that the scenario proposed here can be helpful to think about 
non ergodic dynamics in glassy systems, in particular those in which quenched disorder is 
a priori absent [Q. We note in this respect that the spontaneous appearance of power-law 
distributed 'trapping' times with x crossing 1 at the glass transition has been reported for 
hard-sphere systems [Q. 
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Figure Captions 

Figure 1 Schematic views of the free energy landscape with associated Parisi trees. 1-a: 
REM landscape, 1-b: Full RSB landscape. 

Figure 2 2-a : Fit of TRM decay for Crin at 10 K against two level tree theory with 
simple FDT. Note that the fit is perfect over the whole time domain. 2-b : Fit of TRM 
decay for AgMn at 9 K against two level tree theory with simple FDT. 

Figure 3 Fit of TRM decay for AgMn at 9 K against two level tree theory with gener- 
alized FDT as given by equation (26). 

Figure 4 Fit of the decay of the out of phase susceptibility of CrIn (at T=12 K, = 10^^ 
Hz.) as x"(u;,t) — Xoq. (^) t^^~^, and X2 = .79 (see Table 2). 

Figure 5 Sketch of the temperature cycling protocole and the evolution of x" during 
this cycling (from [|35|). Note the strong spike iust after the temperature decrease followed 
by a perfect 'memory ' in the third stage. 

Figure 6 Aging in n = 1 (approximate result) and n = oo (exact result) ferromagnets. 
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Note that the small time behaviour is regular, at variance with the experimental data on 
spin glasses (Figs 2-a, 2-b). 

Figure 7 M{t^ + t, ty,) log^'^^H) as a function of for AgMn at 8K. The values of 

r = 10~^ sec. and ^ — ^ have been chosen to obtain the best rescahng. 

Figure 8 ^^'^q^^^'^-^^ as a function of for Crin at various temperatures. The values 

of t(T) sec. and qeaIT) have been chosen to obtain the best rescaling of the late part of 
the curve, from which we extract ^ ~ 1.2. 

Tables Extracted values of xi, X2, Qea and ^i, from CrMn and AgMn TRM decay data 
(using simple form of the FDT). 

Appendix A - Some properties of exponential times 

Here we shall review and give proofs of some of the properties of exponential times used 
in this paper. The probability density function of an exponential random variable T of rate 
E, is given by 

p{t) ^ Ee^p{-Et), t>0. 

Property A-1: The exponential time is memoryless, that is conditional on the time T not 
having occurred at some time the distribution of the subsequent time T' before it occurs 
is the same as that for T. The proof uses Bayes' theorem to show that, the conditional 
probability density function for T' is 

Ecxp(E(t + ^ / ^ /X 

and hence the result. 

Property A-2: Given two independent exponential times T and T' of rates E and E' 
respectively, their minimum is distributed as an exponential random time of rate E + E' . 
This can be seen as follows: if the probability density function for min(T, T') is p{u), then 

p{u) = j dtdt'EE' exp{-Et - E't')S {u - min(t, t')) = {E + E') exp (-(E + E')u) . 

By induction we see that for any finite number of independent exponential times the minimun 
is distributed as an exponential with rate given by the sum of their rates. 
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Property A-3: Given two independent exponential times T and T' of rates E and E' 
respectively then 

E 



p(T <T')= f dtdt'EE'e^p{-Et - E't') 
Jt<t' 



lt<t' ' E + E' 

This property generalises to: for j independent exponential times Tj 1 < ? < J of rates 
Ei 1 < i < J respectively, then 

P{T, = min{T„ 1 < ^ < j}) = -3-. 

The proof is simply uses the basic form of property (A-3) and also property (A-2). 
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